Exact solution of generalized Schulz-Shastry type models 
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A class of integrable one-dimensional models presented by Shastry and Schulz is consequently 
extended to the whole class of one-dimensional Hubbard- or XXZ-type models with correlated 
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I. INTRODUCTION 



Integrability has been considered to be one of the most 
striking properties of a model for some time. How- 
ever, the interest was confined to an abstract level since 
integrable quantum systems typically live in one spa- 
tial dimension (but they can. be mapped onto classical 
1 + 1 dimensional systemsEm) . In the meantime, one- 
dimensional systemS|-|have become experimental reality 
(Quantum Hall barsa, polymerso, charge density wave 
systemsu) and thus integrable models have gained imme- 
diate relevance for real physical problems. The Coordi- 
nate Bethe Ansatz (CBA)tljZl solvability is the first mile- 
stone in proving integrability. The (Quantum) Inverse 
Scattering Method together with the algebraic Bethe 
Ansatz complete r-jthe CBA procedure to construct in- 
tegrable modelsM. However, the interrelation between 
solvability by algebraic Bethe ansatz and CBA is still 
controversial. But solvability by algebraic Bethe Ansatz 
is a strong hint for solvability by CBA within an appro- 
priately chosen basis. In the case of the quantum impu- 
rity problems, for instance, the standard CBA procedure 
with plane waves fails, inj-contrast to using a basis con- 
sisting of solitary wavesoEl. The reason for this is that 
the bulk theory already describes fully interacting elec- 
trons and there is no (free) single particle description of 
the system. As a result the scattering with the impurity 
is "diffractive" in the plane-wave basis. Instead, it turns 
out to be factorizable if the many-electron wave function 
is written in terms of kinks and anti-kinks. 
A breakthrough in the theory of integrable models was 
the solution of the Hubbacd Model (HM) obtained by 
E. H. Lieb and F. Y. WuK The HM is a model in- 
cluding on-site Coulomb interaction for electrons moving 
in an atomic lattice. It is believed to,-capture important 
features of high-Tc superconductivitytil and has a metal- 
insulator transition. 

In a recent paper, H. J. Schulz and B. S. Shastry found 
a new class of spivable one-dimensional Hubbard and 
XXZ type modelst3. The modification of the original HM 



and XXZ model consisted in a configuration dependent 
unitary factor in the hopping term. This can be inter- 
preted as an interaction of the charged particles with a 
gauge-field, generated by the density of particles. The 
structure of the unitary factor was expfiiV], where TV is a 
mono-linear functional of particle-number operators; we 
term such models "single particle correlated hopping (1- 
CH)" models. 

The idea behind Schulz' and Shastry's approach is finding 
a basis (through a unitary transformation of the original 
Fock basis) in which the model takes the form of the 
original Hubbard or XXZ model up-to boundary twists 
which do not affect their solvability^. We point out that 
this is equivalent to equipping the plane waves entering 
the CBA with phase factors canceling exactly the con- 
figuration dependent gauge fields in the hopping term. 
In the present paper we will generalize such an idea to 
consider hopping in which N is a multilinear functional 
of particle- number operators. We shall call the resulting 
models n-CH models. We will answer the question which 
n-CH-Hubbard/XXZ models can be mapped unitarily 
onto a corresponding uncorrclatcd but twisted model. 
This finally proves solvability of the model. Conversely, 
a non-removable correlated hopping destroys solvability 
by CBA, since the S matrix becomes configuration de- 
pendent (see Appendix |A|) . 

The paper is laid out as follows. In section [n[ we will 
study the most general form of 1-CH HMs, hence in- 
cluding Shastry-Schulz models. Most of the features of 
the general problem already occur here. Section III ac- 



counts for conserved quantum numbers other than the 
spin orientation of electrons. Multi-chain models where 
the chains interact with each other exclusively gauge- like, 
will be considered. The central results obtained in sec- 
tions H and HI will be used to discuss higher correlated 
hopping in section IV. The 2-CH will be treated in de- 
tail, enlightening the approach to general n-CH. Finally, 
conclusions will be drawn in section |v|. 

A special class of unitary transformations was used to 
remove the gauge-like correlation terms from the hop- 
ping. In Appendix |^ the effect of the complementary 



1 



class of unitary transformations will be elaborated, al- 
ways restricting on automorphic mappings on the class 
of twisted n-CH models of Hubbard- or XXZ-kind. Two 
propositions needed in this section are proven in Ap- 
pendix [5| (see also Ref. |3|) . 



£mT We can further choose £mm = (a non-zero Cm m 
can be included in the parameter Ci,m)- 
We locally transform the Hamiltonian by U: c 



U Cj a U 1 . Number operators remain unchanged, but the 
hopping term is altered 



II. SINGLE-PARTICLE CORRELATED HOPPING 

In this section we discuss a simple generalization of 
Schulz-Shastry models. Such models have the following 
Hubbard-type Hamiltonian 



H 



i 



h.c. 



(1) 



where {cj^cj^,} = S a ^>S^i, {cj >l7 , c;, CT <} = 0, and n ltC , := 

c\ a ci t „. The parameters t, V are the hopping amplitude 
and the Coulomb repulsion respectively. The class of 
models (@) is a generalization of Schulz-Shastry models, 
since a) the parameter A occurs, which means correlation 
between particles with the same spin orientation, and b) 
the spin- and-coordinate dependence of A , a , 7 is unre- 
stricted hereO-i 

We point out that Ajj + i can be set to an arbitrary value 
for all j, without affecting the physics of the model (since 



ni l(T £ {0,1} 



0). We name parame- 



ters like Ajj + i "irrelevant" . A similar argument holds for 
the parameter Ajj: contributions from n JjCT arise only if 
rtj :CT = 1 because of Cj.aiij^ = Cj^ a . Hence, this term can 
be included in the parameter jj (er) . Parameters like Ajj 
will be called "subrelevant" throughout this paper. Irrel- 
evant as well as subrelevant parameters appear as soon 
as phase-correlations among particles having the same 
spin orientation as the hopping particle are involved. 
It is worthwhile noting that Hamiltonian ([l]) is not di- 
agonalizable by direct CBA since the scattering matrix 
is configuration dependent (see Appendix [A|). This de- 
stroys the factorizability of a many-particle S matrix into 
two-particle S matrices. Thus, we first remove the phases 
in the hopping term of ([l]) by a unitary transformation 
and then, we can diagonalize the transformed Hamilto- 
nian by CBA. 

The ansatz for the unitary transformation is achieved 
through the operator 



C j + l,a C 3,<7 



u t 
— * c 1 



j+l,a^3,<? 



exp 



£j",'m) n ™:M 



exp 



i(Cj + l,cr - (j 



2£ 



j,3+V 



The unitary factors in the hopping term of |TJ) are re- 
moved if 



2 (C™ CT ~ £3 + 1,™) = «J,m( Cr ) 



(Cj> - C?+1,<t + = lj(cr), 



(3) 

(4) 



and 



2($£ - gfrj = ; m £ {1, . . . , L} \ {j, j + 1}. 

(5) 

For periodic boundary conditions (PBC), Eqs. (||)-(|^) 
for j = L represent the jump across the boundary. Ad- 
mitting boundary phases in company with the bound- 
ary's crossing, Eqs. (||)-(||) are modified as 

2(C m CT - CD = "i,mW - (6) 



(a,<r-Cl,«r + ^)=7i(ff)-0W- 



(7) 



and 



- ££) = 4» - ^ (*) ; m £ {2, . . . , L - 1}, 



(8) 



In fact, the relations @-(^) and (p)^([|) constitute a sys- 
tem of recursive relations for ^J'^ and Q <cr . 
We will now discuss the exclusions in Eq. (||). The 
corresponding part of the transformed hopping term for 



j is (since £ 



a. a 

3,3 



0) 



! Cj, a exp 



Aj t j(tj))rij !(7 . This term is non-zero only if nj >a — 
1. Hence, this "correlation" factor is equivalent to 



c j+l,<r c j," ex P 
C \+l.<j C 3,° ex P 



For m = j + 1 it is 



[ ( 2 £J-J+1 ~ A 3,3 + l(< 7 )) n 3 + l, 



This term 



U :— exp[i(^ m ni^n mtU + ^i tfl ni tfl )\ — : exp(iS), (2) j s non _ zero on iy if nj+i iCT = 0. Hence, this "correla- 



te use the sum convention) where fy, Ci,-m £ IB arc 
unknown variables which have to be fixed for cancelling 
the unitary prefactor in the hopping term of (|l]). Since 

an antisymmetric part in the parameter £jj vanishes af- Ajj{a) + 2£j£ +1 enters as a modification of Eq. (|f|) for 



tion" factor is equivalent to cj +1 CT Cj.o- irrespective of what 
value Ajj + i(a) — may take. As a consequence, 

there is no condition for m = j + 1 in Eq. (H), and 



ter summation, it can be defined fully symmetric: £ 



7j (see Appendix |B|) . 
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FIG. 1. Elementary viable (leading to (|ic|)) and 
not elementarily viable (called "exceptional" ) loops, are 
drawn in the space of indices of (the crossings 
of the dotted lines are the possible index points). The 
expceptional loop determines the boundary phase for 
the parameter A. Full arrows mean an application of 
the recursive relation (^), which can transport in hor- 
izontal direction only. The other arrows visualize the 
application of the symmetry £ij = £j t i. The double 
lines are reminders of the missing connection between 

\ j an d Cjj 7 - Hence, it must not be crossed by the 
full arrows. 

The system of equations cannot be solved for 

arbitrary a and A (the number of parameters on the 
right-hand side is larger than on the left-hand side). To- 
gether with the symmetry of Eqs. (||) and (g) de- 
fine the effect of an increase in both site indices of ^.j; 
namely an increase of £ by the related a respectively A. 
Starting from an initial parameter, say £i i := 0, every 
£ij is then defined by passing from (1,1) to and 
summing up the contributions from the recursive rela- 
tion. £jj is well defined if this procedure is path inde- 
pendent. This is equivalent to demanding that contribu- 
tions from closed loops in the plane vanish. To ver- 
ify this, it is sufficient facing the smallest possible loops: 
(j,m) -> (j + l,m) -> (j+l.m+l) -> (j,m+l) -> (j,m) 
(which we will henceforth call "elementary"). Applying 
Eq. (|) respectively Eq. (|), we obtain 



ocj, m +i{o) - atj, m (cr) = a mi j + x(—a) 
and 



,'(-0, (9) 



A hm+1 (a) - Aj >m {&) = A m j + i(a) - A m j(o) (10) 

for m ^ j,j±l. We call the conditions (||), (|l(])) "closed- 
ness conditions" . The recursive relation and the closed- 
ness condition can be written in a more compact and 



clear form in terms of the discrete gradient, defined by 
d x f(x) := f(x + 1) — f(x). The recursive relations (g) 
and (H) then read 

- 2d x t?'- a {x, y) = -2d x C a ' a (y, x) = a(x, y; a) (11) 
-2d x ^(x, y) = -2d x ^(y, x) = A(x, y; a) (12) 

and the closedness conditions take the form 

d y a(x,y;a) = d x a(y,x;-a) (13) 
d y A(x, y\ a) = d x A(y, x; a) ; x ^ y. (14) 

With these conditions being fulfilled!!!, the correlations 
from the hopping term can be removed and the rotated 
model is finally known to be solvable by CBA. 
For open boundary conditions the correlated hopping can 
be "gauged away" completely, yielding the HM without 
any boundary phases. 

Instead, PBC lead to the HM with twisted boundary con- 
ditions. Periodicity implies that the parameters a, A 
and 7 are periodic in their site- indices with period L. The 
boundary phase is determined by hopping from site L to 
site L + 1 = 1; such a phase enters as written in Eqs. @- 
(|) for the relations £i, m <-> £i >m and Ci,m <-» Ci,m and 
then for A\ m <— > AL mi Oi\^ m <— > ctL.m- Equations (§)-(§) 
lead to 

L-1 

2 (£i,m CT _ S; m ") = E a h™( a ) = 

3=1 



-a L , m (a)+(t>^{a), (15) 

3=1 

1 VSm-l,m 5m,m+l/ 
^-Ai.rnH+^ftV), (16) 



L-1 



L-1 



Ci,<7 - Cl, ct = E( 7 J'( CT ) + A i'^ 

3=1 

= -Tl(o-) - Al,l + </>(o-), 



(17) 



where <fi denotes the boundary phases. They can be de- 
termined (without solving for the £'s explicitly) from Eqs. 
Q to (01 as 



4? 



a(!) 



( CT ) = E a J> m (°')' 

L 



(18) 



|| i 17 ) ~ E A 3^( a ) + 

3=1 

l(f)+An-l,in+l( (r )i ( 19 ) 

L 

<^)=E(^( CT ) + (2°) 

i=i 



3 



Eqs. and ( |To| ) also ensure that the phases in Eqs. JTs| ) 
and ( |19[ ) are TO-independentO. Summarizing, iff Eqs. 
(^) and ([n]) are fulfilled, the Hamiltonian ([j]) can be 
mapped by U (Eq. @) onto the usual HM with modi- 
fied boundary conditions. The boundary twists are given 
by 



:= <P(<j) + <f>f}{a)N- a + $h<T){N a - 1). 



(21) 



We emphasize that here, the factor (N a — 1) appears in- 
stead of N a . The reason is that the recursive relation (||) 
does not exist for m = j,j + 1. But since one particle 
with spin a has to be on site j and none on site j + 1 (or 
vice versa: see the corresponding hopping term in ([[ 
one particle less accounts for the phase. 



A. Translational invariant models 

In this section we assume a translational invariant 
model, and hence translational invariance of the param- 



eters a and A: a 7 - Aa) 



-,(ff), A u (a) =A j -i{a). 



Restricting ourselves to translational invariant unitary 
transformations!!^ [£,j t k — £,j-k), the former elementary 
loops become multiple loops (see Fig. |J). Then, the ele- 
mentary loops are (j, m) — ► (j + 1, to) — * (J + 1, to + 1) 
and the closedness condition only consists of the follow- 
ing terms from Eqs. ([))) and ( |l0| ) ctj(a) = — ct-u+i) (— cr) 
for all j, Aj(er) = — A_y+i)(o") for all j ^ 0. Also here, 
for the parameter A not all elementary loops are viable. 
But the only exceptional loop emerges from the case cor- 
responding to to = j in Eq. ( ^0[ ) (that is a square on the 
diagonal - the fat dotted line in Fig. Q). But as already 
mentioned, the (excluded) condition for to = j in Eq. 
( |io| ) is trivially fulfilled. This is confirmed by calculat- 
ing the effect of the jump from to 2 in the exceptional 
loop in Fig. [I], which here constitutes an elementary loop 
(comparing with Fig. ^ it turns out that the points and 
2 are identified since translational invariant parameters 
are assumed). As a consequence, the feature of addi- 
tional phases imported by subrelevant parameters and 
hence non-viable loops are absent here. 
But we want to stress that even here Eqs. (@), (]l(]) are 
necessary and sufficient conditions for solvability - for pe- 
riodic as well as for open chains. This is because a trans- 
lational invariant model need not necessarily be gauged 
away by a unitary transformation with translational in- 
variant parameters. 

Now, we can connect to Ref. |l2|. We obtain the bound- 
ary phases in Ref. [l^ through the identifications ot^ m ~ > 



otj-m and otj(a) 



-craj 



In addition, A, 



j,m(o") - 

7j(er) = and flij — a/2 ■ giving an antisym- 

metric Pij. However, in order to guarantee closedness 
and hence solvability of the recursive relations, 



OLj— fji Otj—fji—\ 



-(0t r , 



(Xm-i-l) 



has to hold. 



A special type of such models is studied in Ref. [l5| 
with the only non-zero values «o = a -i = Vi where rj is 
a real parameter. The closedness conditions are all ful- 
filled and therefore the correlation factors in the hopping 
can be gauged away giving exactly the phases found in 
Ref. However, we could not reproduce the results 
following the path suggested in Ref. |lj (see the general 
discussion in appendix |A|) . 



III. MULTI-CHAIN MODELS 

Multi-chain models can be treated analogously. The 
chain-variable plays the role of the spin-variable now. 
The only difference is that in general the chain variable 
ranges in a set different from { — 1/2, +1/2}. Since the 
method we employ to diagonalize the Hamiltonian stays 
the same, we only state the results. For multi-chain mod- 
els we assume intra-chain hopping only. 

At first, we deal with spin- less fermions and let the 
hopping-term of (Q) take the form 



a c,j+i a c,j ex P(i7j(c)) exp[L4^W ; ] + h.c. 



(22) 



where c and d are chain indices, which run in {1, . . . , C'} 
with C being the number of chainstj. The unitary trans- 
formation will be 



U := expfi^m^m^™ + CiM»fy,i)] , 



(23) 



where fj,, v are now site indices (instead of spin indices of 
Sec. II), running in {1, ... , C} (instead of {-1/2, 1/2}). 
In this case, the recursive relations are 



2(^f m - dt h J = A^ m ; to ± 3,3 + 1, (25) 
(Cy - Cc J+ l + = 7cj , (26) 

The closedness conditions read 



" A fm = ~ ; (C ? d) V TO + j ± 1, (27) 



a; 



-A 



; c ^ d V j ? 0, -1 



for all c, g?, j, m respecting the noted exceptions. Equa- 
tion ( p8| ) applies to the translational invariant ansatz. 
The boundary phases finally become 



* c = 0c + ^^ 1) (d)(A^-<M> 



(29) 



where 



(d) = J2 A 7^ c ^ d ' 

3=1 



(30) 



(28) 



4 



L 

E 



A" 



A 



m— l,m+l ' 



L 

E 



A c - 



(31) ^W(c,<7) 



(32) 



To compare this with Rcf. one has to replace 
€l,m witn k^vB,j,^{l - to) and additionally A£ m with 
— £,c,dA Ct d(j — m). There, £ Mj „ was antisymmetric and 
hence, B^^il — to) has to be antisymmetric, too: 
B^vim) = —B Vjfi (—m). The solvability conditions ( ps| ) 
then transport into 

Cc,d^. c ,d(j - m) = -£d, c Ai,c(ra - j - 1) 

^c,d(j - m) = A diC (m 

for all to, j, c, d. This is equivalent to A c ,d(— m) = 
Ad, c (jn ~ 1) for all to, c, d as obtained in Ref. EOl This 
condition is sufficient for CBA solvability but not neces- 
sary Necessary and sufficient is A c _d{— to) — Ad^ c (—m — 
1) — A c ,d.(m — 1) + Ad, c {m) — for all to. 

If spin has to be included, nothing changes but the 
number of indices. The hopping term becomes 



exp(i7j (c)) exp [iAj' t (a, p)n d} i- p ] + h.c 



(33) 



L 

E 

3=1 
trn — 1, 

+A' 



.4' 



c,c 
m.m- 



^ c _ l!m+1 (a,a), (40) 



^(7^(«T)+^;J(<7 )( 7)) 



(41) 




where c, d are the chain indices again and cr, p are spin 
indices. The unitary transformation takes the form 

U := exp[i(^(p,T)n Mi/;p n„ :m;T + C,ix,l; P n ^,l;p)] > ( 34 ) 

where /i, ^ are the chain indices and p, r are spin indices. 
The recursion relations read 

2(^,1 w, p) - p)) = ^>>> ^ ( 35 ) 

c / e? V a ^ p V m / j,j + 1, 

(Ccj(CT) - Ccj+lW + 2£^ +1 (cr, cr)) = 7 CJ (cr), (36) 

and the closedness conditions are 



A c /i +1 (a,p) - A^a, p) = A^f j+1 (p, a) - A d m c Ap,a)- (37) 



FIG. 2. On the left, an elementary loop in 
two-coordinate space is drawn. The circles are the co- 
ordinate sites, say (i, j), corresponding to a parameter 
. If the parameters in the transformation are as- 
sumed depending on the difference of the coordinates 
only, the dotted lines connect identified sites in index 
space, corresponding to identical parameter values. One 
can imagine the resulting reduced lattice by furling the 
coordinate space along the dash-dotted line to a cylin- 
der or by pushing the identified sites on each other, 
which is visualized on the right. In this case, the for- 
mer elementary loop turns into a double loop. Each of 
these two loops is an elementary loop. 

The overlap of Rcf. |lj with the work presented in the 
Refs. p2p^ ,mas 1 discussed by the authors and was argued 
to not existEJo. Now it is shown that the correlation in 
Refs. ^2[^3| cannot be gauged away, since the correlation 
parameters violate Eq. (|37|). However, the Hamiltonian 
having the slightly but essentially different phases 



exp 



ni— 1 - 



2tt 1<7 1 (42) 



N a 



c^d V cr^p V m^j±l, 

for all c, d, j, to, cr, p respecting the exceptionsill. 
boundary phases become 



The 



instead of 
exp 



n l+l, 



ni-l. m ,o \ + 



2^7^ 

Nn 



} (43) 



where 



c 

E E 

d=ipe{U} 



yields the BE found in Reference E 
The models studied in Ref. 



are equivalent to multi- 

^lid, p)(Nd. p — Sc^Sp.a), (38) chain models of spin-less fermions. In the first model, 

, , . tin m — 3— 1 



the only non-zero model parameters are A 



-A 



m+l,m 



-46 and 7j(to) 



m,m+l 
3,3 

2(6m+l .in 



6« 



(d, P )=j2 A 7.U°>py> c + d v ° + 



(39) 



l 3,3+l 

©Tn,m-i)j where we used the notation in Rcf. |2d: 
j and to are the site and chain index respectively. 
The second model is represented by the parameters 
• 1 y , . 1 =-A jj = B m , m+ i+a m , m+ i and A jd = 



5 



-A 



in , m + 1 



e, 



tions are fulfilled iff O m ,m+i 



-i. The closedness condi- 
= ©m+l,m and am,m+l 



IV. TWO-PARTICLE AND HIGHER 
CORRELATED HOPPING 

The procedure developed in the previous section can 
be extended to consider higher (than one-particle) corre- 
lated hopping. First, we face explicitly 2-particle corre- 
lated hopping (2-CH). Then we will sketch how to deal 
with the general case of n-CH. 

2-CH corresponds to the occurrence of a term like 
u,'it{j,cr)ni l \n m , l j, in the exponential factor of the hop- 



ping term of (|l|) 
is 



The 2-CH Hubbard-type Hamiltonian 



H = 



3,o 

x exp 



x exp [ia^(j, <j)n m ^\ + h.c. 



(44) 



Without loss of generality the parameters d^(j, cr) can 
be chosen symmetric in the index pairs (I, A) and (to, a) 
and vanishing if these index pairs coincide (see Ref. 26 ) . 
The parameters cij+i m (j, cr) and aj +1 (j, a) are irrelevant 
for all j, to, cr, fi; the effect of the subrelevant parameters 
on the lower correlated ones will be discussed later on in 
the present section. 

We first remove the phases in the hopping term of 
( [f~4| ) by a unitary transformation. Then, we diagonal- 
ize the transformed Hamiltonian by CBA in computing 
the boundary phases. The 2-CH demands an exponent 
^m'r ni ^ nm ^ nn -p m ^ ne unitary transformation U : 



U 



exp[i(f* p 



(45) 



Both £ and £ are totally symmetric and vanish if any two 

pairs of parameters coincide. 

The hopping term is transformed into 



t u t 

c j+l,o c 3,o > C j + l,<r c j,<? X 

x exp[3i(Cj; A i^ m - Q^)ni :X n m ^ + 

whereas the Coulomb interaction term remains un- 
changed. This leads to the recursive relations (compare 
with (§- ©) 



Off cr,A,/i 



) = (46) 



(to ^ Z V /i ^ A) and (Z, m ^ j, j + 1 V X, fi ^ a) 



2(£ 



(T,fl 

j,m 



(47) 
(48) 



£7+1,171 "I" ^^jj'+l,m) — a m(jj CT ) 
m ^ J: j + 1 V /! ^ a 

(Cj> - G+i,<7 + 2 £j°j+i) = 7i(o") 
for the parameters in /7. 

We point out that, in the present case, two kinds of el- 
ementary loops exist because of the variety of indices 
in ^m0»- Namely: (j;l,m) -> (j + l;Z,m) -> 
(j + 1; Z + 1, to) -> (j; Z + 1, m) -> (j; Z, to) and (j; Z, m) -> 
(j + l;Z,m) -> (j' + 1;Z,to + 1) -> (j;Z,m + l) -> (j;l,m). 
However, due to the symmetry of the d both loops give 
the same closedness condition for a, which is 



a 



A . fx 

l+l,m 



-CT.LL 
V 



(Z,A) 



a 



(Z,A) 



(49) 



for Z ^ j, j±l V A^d and j,j + 1(1,1 + 1) V ^ ^ A. 
The corresponding boundary phases are 



L 



(50) 



3 = 1 



^ cr, 
L 



(51) 



3 = 1 
jVM±1 



+ai-!,i(l + l,cr) + a^^C/.a) + 



(52) 



where Z 7^ m, to ± 1 can be chosen arbitrarily. The result 
turns out to be independent of this choice. For less than 
three sitescJ, Eq. (p2h is ill-defined. It reflects a physical 
limitation: for the boundary phase ^ (cr, cr) to occur, at 
least three particles with the same spin orientation have 
to exist; this is possible only if at least three sites are 
available. 

Now we will discuss the effect of the subrelevant param- 
eters a.j'£(j, cr) and aj(j, cr). The 2-CH subrelevant part 
of the exponent in the hopping term is 2(dJ'^(j, <t) + 
S^j'j+i m )- As discussed in the previous section, this term 
does not vanish in general!! because the recursive rela- 
tions do not cover the index grid of £ completely. It 
contributes to the 1-CH part instead. It has to be added 
to the right-hand side of Eq. (|47|). As a consequence, the 
parameter £ drops out and the recursive relation reads 



(53) 



G 



where /3^(j, a) is defined as 



(54) 



Doing the same for the 1-CH subrelevant part in the hop- 
ping (concerning Eqs. (^) and (^3|) now), C, drops out 
and in Eq. ( f48| ) , (cr) will be substituted by (<r) 

= 7iW+<0». (55) 



Using Eq. ([53|), the second set of closedness conditions 
are obtained 

fa+iti, *) - v) = fi) - P°(m, M ) (56) 

and the boundary phases are 

= 4><r + ]T (A)(iV A - S X ,a) + 
A 

+ ]T<^ 2) (A,u)(jV A - SxJiNp-Sw), (57) 



where 



3=1 

L 

W(a)= £ + 

3=1 
3 ym-l,m 

/?^_ 1 (m,<7)+^ +1 (m-l )( 7), (59) 

(60) 



i=i 



The analysis presented above can be generalized to con- 
sider n-CH. 

(n) (™) 

Let a , (3 be the parameters of the n-CH part of the 
hopping: 



(n) 



/3 {i}0» := 



n + fc\ (n+fe) 



fe=0 



(61) 



a <*) + (« + i)^ 1 ^!} 0'i ( 62 ) 



(n)° 



(n)° 



/?,,... (J» := /?,,,.... (.,0 :=0. 



(63) 



(for multi-chain models the spin index is a multi index 
spin/chain). The dots (...) stand for an arbitrary series 
of indices except they appear in between two equal in- 
dices (a . . .a). In this case, the void is meant to be filled 
up with k times a. {S} is a set of spin- or color indices, 
{1} a set of coordinate indices. The sum in Eq. ([H]) of 
course is finite. The variables are confined since CL — 1 
is the highest possible correlation level if C is the number 
of inner degrees of freedom (C = 2 for spin 1/2). Hence, 



k < CL — 1 — n and for n-CH with n < CL — 1, the sum 
in Eq. (|6l| ) already stops at k — n — n. Furthermore, 
as mentioned above, parameters with coinciding index 
pairs can be assumed to be zero. This leads to Eq. (|6^). 
Eq. ( |63| ) only reminds us that all subrelevant parts are 
removed if (3 is used instead of a. 

The closedness conditions take the form 



i+i, 



The n-CH boundary phase is given by 

= $C"-i) + J- 4»)(A, . . .)(N X - 8 Xt(T ) 

A,... 

where the phases <p n > are 



(64) 



(65) 



(n) 



L („){«} 

({£}) = £ (3 {J} (?»; 

3=1 

a £ {5}, 

L (n) <7,...,<7,{S'} 



(66) 



3 = 1 
i=l..k 



A r (n) 



(n) c 

/3. 



...,<r,{S'} 

,1*— 1,. ..,{!'} ih,cr) + 
r,...,a,{S'} 

<7 {5'}. 



(67) 



The pair of index sets {5}, |{<S}| = n and {X}, {I}|=n, 
respectively {5'}, |{<S'}| = n — k and {I'}, \{I'}\ = n — k 
must not have coinciding index pairs and in Eq. (67) no 
two U must be identical or neighbored. Thus, Eq. ( |6"7j ) 
holds for L > 2k. This validity range can be maximally 
enlarged using the closedness conditions (see Eq. (0)) 



(n)° 



.,<x,{S'} 



b^(a,...,a,{S'})= J2 0l+l,...,l+k,mti,<r) + 



3=1 



[/,/■ 



fc (n) <T,...,<T,{5'} 

i=0 

k + l] \ {l + i,l + i + l}, 
a $ {S>}. 



(68) 



This formula holds for L > fc + 1, which is a limit 
set by physics - analogous to the 2-CH case. The re- 
sult is ^-independent. We point out that exp[LF({n})] 
with F({n}) being an arbitrary functional of number 
operators is not the most general unitary operator in 
Fock space. The most general is exp[iG({V, c})] where 
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G({c^, c}) constitutes a Hermitean functional of the com- 
plete set of creation/annihilation operators. In the Ap- 
pendix |c| we will show that the class of unitary opera- 
tors discussed so far is large enough as far as removals of 
phases in a nearest neighbor hopping term are concerned. 



V. SUMMARY AND CONCLUSIONS 

In summary, we found a complete characterization 
of Coordinate Bethe Ansatz (CBA) solvable Hubbard- 
type Hamiltonians with unitary correlated hopping for 
fermions. A necessary and sufficient criterion for such a 
Hamiltonian being solvable by CBA was formulated (see 
Eqs. (|), ©, and ©). 

In contrast to what is suggested in Ref. [15], we find 
that these models are not CBA solvable in the ordinary 
plane waves basis. Indeed, in such a basis the scatter- 
ing matrix is configuration dependent (see Appendix ^) 
thus describing diffractive scattering. The particles inter- 
act non-trivially even for vanishing Coulomb interaction. 
Such a situation resembles the case of the impurity prob- 
lem in the sense that also there, the free picture already 
contains some residual interaction due to the impurity. 
Solvability of the models is recovered if the correlations 
from the hopping terms can be gauged away. This is 
equivalent to equipping the plane-wave basis with ad- 
ditional density dependent phases. Only in this modi- 
fied basis can the correlated hopping be absorbed in a 
boundary term. For the models considered in Refs. |2^]23| 
no such basis exists and hence, they are not solvable by 
Bethe ansatz as they stand. They can however be re- 
paired by modifying slightly their hopping term, as done 



Refs. [13L[29 30. Even for absent Coulomb interaction, 



at the end of chapter III 



The boundary twists for solvable models with periodic 
boundary conditions are given explicitly in this work. 
The corresponding Bethe equations are known from Ref. 
|l3| adapting the boundary phases for a spin degree of 
freedom only 



H ifapj -co + m 



e ip * L = c 



fjr i((a -Cb) + Tt _ p -i(* T -*i) rr l( - sinpi ~ 



V 
21 



n i( S mpi ~ Ca) 



(69) 



41 



41 



Therein, pj are the quasi-moment a in the plane waves 
used in the Bethe ansatz (see Eq. (Al)), Ca are spin ra- 
pidities, t and V are the Hubbard model parameters and 
$ CT are the boundary phases, which have been determined 
in this paper. The connection to Ref. |l3| shows that in 
the solvable cases, the correlation in the hopping term 
of the Hamiltonian (see, for instance, ([[])) is equivalent 
to applying a magnetic flux to the system. However, 
such a flux is generated by the particles themselves (in 
particular, it is not an external magnetic flux). Ground 
state properties can be deduced from those calculated in 



the many particle energy is not a sum of single particle 
energies. The effect in the energy density is of first order 
in l/L and in — $jj in the thermodynamic limit. So 
one can argue that correlated hopping accounts for a non 
trivial interaction between the particles even for vanish- 
ing Coulomb interaction. 

The results obtained here can be applied to models for 
particles with deformed exchange statistics(DES). This is 
done via a mapping from DES to CH modelsLj, whereas 
special DES modelsphpje been discussed recently using 
direct Bethe ansatztaffl The details will appear in a 
forthcoming paper. 
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APPENDIX A: CONFIGURATION DEPENDENT 
S MATRIX FROM STANDARD CBA IN 
PRESENCE OF CORRELATED HOPPING 



In this appendix we explain why correlated hopping de- 
stroys solvability by direct CBA for Hubbard- and XXZ- 
type models. 

The CBA is an ansatz fop-.the wave function in a so- 
called fundamental regionO, in which the interaction 
term doesn't contribute. This fundamental region has 
to exist at first. For the Hubbard model, this is guaran- 
teed by particle-hole symmetries. This gives the energy 
in terms of the distinct quasi-momenta. The wave func- 
tion has to be defined uniquely on the intersection lines of 
the fundamental regions, where for the Hubbard model 
also the interaction enters. Both demands yield the S 
matrix, which represents the effect of the interaction, a 
scattering of two particles. For Bethe ansatz solvability, 
the interaction must not have an effect beyond permuting 
the quasi-momenta of the particles. This demands that 
the S matrix fulfills the Yang-Baxter equation. Now the 
boundary conditions, if compatible with the S matrix, 
lead to additional conditions, fixing the eigenfunctions 
constructed by Bethe ansatz up to normalization. As 
a consequence, any additional condition destroys Bethe 
ansatz solvability. Every correlation in the hopping fur- 
ther restricts the parameters in the Bethe ansatz. 

Let us assume a given CH Hubbard-type model can be 
transformed iso-spectrally to another CH Hubbard-type 
model, for which fundamental regions exist. This should 
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mean that neither interaction nor correlated hopping con- 
tributes. All correlations already contribute, when no in- 
teraction is yet to be included. The Bethe Ansatz will 
have the form 

N 

ip( Xl ,...,x N ) =: 2^ ^(7r|7r')e »=* (Al) 

where the permutation 7r' is chosen such that an appro- 
priate order in (xi, . . . , xjv) is achieved. Here, the cho- 
sen order which defines the fundamental regions, will be: 

X\ < . . . < X ]\f . 

1. "Pinned" correlations 

At first assume a correlation which appears if a par- 
ticle sits at a special site. Let particle number j sit on 
this site. The effect of this is a shift in all momenta 
by the correlation strength ip, except that momentum of 
the designated particle. This is seen by projecting the 
Schrodinger equation on the specified configuration. The 
particle causing the correlation with the others feels no 
correlation from the hopping term applied to it, since it 
will be transported to that site by the hopping. Thus the 
resulting term for the energy is 

N 

E = -2t^cos(p 7r(i) +ip)- 2tcos(p 7r(j) ). (A2) 

In this equation, 7r is the momentum permutation from 
the Bethe ansatz. This energy is neither independent of 
the permutation as it ought to be, nor does it coincide 
with the original energy formula E — — 2t cos(pj). 
Note that ip could even depend on the spin orientation 
of the considered particles. 

2. Relative correlations 

Assuming the hopping term to be 

c ]+Im c j,<7 exp(i^ra j+A , M ) + h.c, (A3) 

where a and /j, are spin indices. Further assume one single 
particle at coordinate j to be affected by this correlation. 
The corresponding condition from the Schrodinger equa- 
tion is 

^(j-l,a;j', A1 )(e i ^-^ -I) 

+ip(j + 1, a; f, n) (e^+^ - l) = 0, { 1 

where we omitted the spin index from the argument of 
the wavefunction. Two cases can appear independently 
from each other: 

• a particle with spin orientation /i sits at site j' = 
3-1 + A 



• a particle with spin orientation /i sits at site j' = 
j + A 

leading to — 1, cr; j — 1 + A, ^) = in the first case and 
+ 1, cr; j + A, /i) — in the second case. These con- 
straints appear in addition to the usual "continuity" con- 
dition arising from tp(j, a; j, p) = cr; j, fi). If cr = ^t, 
at the most two of these three conditions can coincide if 
A = or A = 1, which however constitute a subrelevant 
and irrelevant correlation, respectively. For /i = — cr, one 
of the correlation terms then coincides with the on-site 
Coulomb interaction. The other, however, is still remain- 
ing. Tlie only way out is to modify the Bethe ansatz 
slightlyEa. In Ref. B3, the fundamental regions are with- 
out intersection for different kinds of particles (different 
spin orientations). Consequently, for different kinds of 
particles the "continuity" condition is absent and hence 
substituted by one equation coming e.g. from correlated 
hopping. The second equation has then to coincide with 
the interaction. This indicates that the hopping term 

C j+l,a C j,a ex P( i( P n J+^ -°) + h - C - 5 A ^ {0, 1} (A5) 

is the only possible relevant correlated hopping which is 
tolerable by direct Bethe ansatz for Hubbard-type Hamil- 
tonians (for XXZ-type models no relevant CH is treatable 
by direct Bethe ansatz). But then, the obtained S ma- 
trix still had to fulfill the Yang-Baxter equation. In Ref. 

, this kind of hopping term is studied. There, ip is 
purely imaginary and A(cr) :— (a + I)/2 G {0,1} was 
chosen. This Hamiltonian is shown to be solvable by a 
modified Bethe ansatz in the absence of Coulomb inter- 
action, V = 0. For V ^ 0, the resulting S matrix no 
longer fulfills the Yang-Baxter equation. 

To summarize, except for a very special type of relative 
correlation, namely with particles at distance zero or one 
towards raising site number, one obtains independent ad- 
ditional sets of conditions which the wave function had to 
fulfill. This obstructs the direct Bethe ansatz, since the 
S matrix becomes configuration dependent. It then no 
longer factorizes into two-particle S matrices. Regard- 
ing the exceptional CH mentioned above, a configuration 
independent S matrix is obtained by applying a variant 
of Bethe ansatzE3. But for finite Coulomb-interaction 
strength V, this 5* matrix does not fulfill the Yang-Baxter 
equation. As a consequence, no CH XXZ-type model and 
no interacting CH Hubbard-type model is tractable by 
direct Bethe ansatz for finite V. 



APPENDIX B: A NOTE ON SUBRELEVANT 
PARAMETERS 

In the following, we will discuss the contribution of the 
subrelevant parameters. Recalling the definitions, irrel- 
evant parameters like Aj_ij(cr) do not at all affect the 
physics of the model, whereas subrelevant parameters like 
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Ajj (a) contribute to the uncorrelated part of the hop- 
ping. For this reason no recursive relations come from 
them. It is worthwhile noting that Eq. ( ]l0| ) is trivially 
fulfilled for m — j. For this case, the exponent produced 
after transformation by U is m m (2£j_^ a +Ajj(a)). This 
term hence appears as a phase additional to 7j(c). The 
parameters Cj+i j are ( U P t° an a dditive constant) given 
by jumping from to 2 in the exceptional loop in figure |l| 



£fij) = ViW + Ai-ij+ii*). (Bl) 



We will now consider the special case in which ^Cj+i j + 
Aj,j(a) = holds. This is exactly what the relation (^) 
would result in for m = j. It implies that the subrele- 
vant parameters do not create additional phases in the 
uncorrelated part of the hopping. Together with (|Bl| ) it 
gives 

AjA*) ~ A i,i-i{°) - ^i-y+iW - Aj-ij-xi?) = 0. 

(B2) 



The condition for the relations (@) for 



J 



and 



= j + 1 being consistent is A 



= -A 



which 



can always be fulfilled by properly choosing the irrele- 
vant parameter Aj—ij. Inserting this into Eq. (B2), one 
finally also exactly obtains relation (||) for m = j + 1. 
We can therefore conclude that if Eq. ([B2| ) holds, the 
contributions from subrelevant parameters cancel out. It 
bridges the void in the recursive relations (^). 
Since all the a's are relevant, no voids occur in their re- 
cursive relations, which is equivalent to all elementary 
loops being viable. 



APPENDIX C: UNITARY TRANSFORMATIONS 
IN THE FERMION U(R) ALGEBRA 



In this appendix, the transformed number operators 
c^Cj and hopping operators c|c^ will be studied for the 
most simple unitary transformation such as 



U = exp 



ia(c' k ci + cjck) ; k^l 



(CI) 



The r 2 operators {c\ci , 1 < i, j < r} span the Lie alge- 
bra u(r): 



(C2) 



where the Cartan basis is generated by Hi — ejej, 
(i = l...r): 

[H l ,c s J c k ] = (Sij - 6i,k)c}c k (C3) 
The u(r) algebraic structure allows us to write 



e k - l i 



- sin 2 a(5i. m - 5 k , m )(ni - n k )+ 

-M(sinacosa - a){8i <m - S k , m )^ki 



(C4) 



fe,m i 



+4 sina[ < i>^ ; ((5fe jm (5; !ri + Si tm S ktn )+ 
where we applied the following equivalence for adjoints 

oo 

[Ad(exp(A))]B := e A Be~ A = ^ —[A, B] n =: exp [ad(A)]B, 

n=l U ' 

[A,B] Q :=B, [A,B] n+1 := [A,[A,B] n ], 



and we have defined <& kl := 



cjc k , and $ fc , 



c\c k (the $ + are Hermitean, whereas the $ are 



4 c i 

anti Hermitean). 

^From this it is seen that the transformation of the Hub- 
bard Hamiltonian creates arbitrary-range hopping from 
and to the sites k and I, as well as pair-hopping cre- 
ated from the interaction term. This can be understood 
from interpreting $? kl as a Hamiltonian itself and a as 
the time. Then, Eq. (C4) gives the number and hop- 
ping operators in the Heisenberg picture. ^From this 
interpretation it seems reasonable that no linear combi- 
nation of <E>^" ; with k I will ever be able to just remove 
phases in a nearest-neighbor hopping term. However, 
this is not absolutely true. To point out the exceptional 
cases, the investigation has to be completed. Since the 
identity [Ad(exp(Aj)]B = exp [ad{A)]B does not consid- 
erably simplify calculating the action of a more general 
unitary transformation, another approach will be taken. 
But at first, the problem will be reduced as far as possi- 
ble. 

In a general product of creation and annihilation op- 
erators, one can at first collect operators occurring in 
pairs using the exchange algebra. The result is a mul- 
tilinear form of number operators besides a multi-linear 
form of creation and annihilation operators, so that no 
two operators have coinciding indices. So the most gen- 
eral exponent appearing in U can be written as follows 



Hk,l;{p},{<l} n k n l Hi c li C H + ■ ' ' • 



(C5) 



Here, the different £ are symmetric in the indices before 
the semicolon and Hermitean in the index sets behind 
it. They vanish if any two indices coincide. Assuming 
a pure m-linear form of the creation/annihilation opera- 
tors, then a transformed bilinear object contains m-linear 
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and even higher terms. This leads to multi-particle hop- 
ping and interaction terms including more than two num- 
ber operatorscS. Since the aim is to stay in the class of 
Hubbard- or XXZ-type Hamiltonians, we consequently 
can limit ourselves to general multilinear forms of num- 
ber operators only, as already studied above, or bilinear 
forms of creation/annihilation operators only. We will 
discuss the latter in the following. 

U := exp[i&,,4c,] ; = 6% , = 0. (C6) 

We point out that t/c^Cnt/ -1 is bilinear, since 



4 



[c' k ci, c\. n c n \ is bilinear (see C2). So, one can determine 
the result by projecting on the desired initial and final 
states. 



, -in 

(°\ c m, p,'4 c ij c m< 1°) = 

= <o| c m/ kkA c i] 4,-! 1°) <°l c 



• ••4, |o> (o|c mi 

sm f ,m n -iCm n -i ,m n -2 



|0) = 



£,m 1 ,m i — (£, )mf, 



all written in sum convention. One directly obtains from 
this 

(0|c rof C/4j0) =(expi£) w (C7) 
(0| c^f/c^t/t^ |0> = (expiO m/ , fc (expiO* ma , (C8) 

where hermitecity of the £ was used. The indices 
fc, ^ are not summed over. It can be shown that the 
transformed interactioxuterm can never include a single- 
particle hopping termE3. This already proves that the 
reverse direction is also impossible. Thus, number opera- 
tors and nearest-neighbor hopping have to remain "type- 
invariant" under the transformation, since the type of 
Hamiltonian should be preserved. This results in restric- 
tions on the matrix £ 

( ex Pi£) m/ ,fc (expi£)m i>fc = V m f + m i ( C9 ) 
(expi£) m/J+ i (expi^^ j = for \m f - rm\ > 1. (CIO) 

The first conditions emerge from transforming the in- 
teraction term, whereas the second one comes from the 
hopping term. Using both, we can deduce the structure 
of exp i£. 



Vfc 3 1 m(k) 3 (expiO m(fc)ifc ^0 
m(k) = k + const. 



(expi£) fei (a) =: 8i, k+A r k e^ ; Y[ k ne^ = 1. 




Unitarity of the matrix implies 



1 for all I. 



Hence it can be assumed that r; = 1 for all I. So, 
one finally concludes that U(a) transforms ct+l c j 

Number operators remain un- 



J 



-i+A+x°i+A- ■ -■— — -x- — 

changed. These phases can be gauged away leaving no 



K<t>3 



boundary phase. It is worth noting that pure d-range 
hopping on an L-site chain can also be obtained iff d and 
L are relatively prime. With the analysis used here it can 
finally be shown that unitary transformations with multi- 
nomials of odd degree always produce particle-number vi- 
olating terms. Hence, they also make us leave the class of 
models we consider. But it is clear from this, that a huge 
class of models, which is far from being Hubbard-type 
can be constructed by unitarily transforming the Hub- 
bard model. They all are solvable and have the same 
spectrum as the Hubbard model. 



APPENDIX D: PROOFS 



1. Proof that transformations such as (C5) exceed 
the considered class of models 

Consider unitary transformations U given by 
U = exp i(Q + 5) 

3 = £{pm?} il 4i c n + &;{p}.{5}"* rii 4,%+ 

Hk,l;{Ph{l} n k n l IL 4i C H + • • • ' 

where is a multilinear form of number operators only 
or bilinear in creation/annihilation operators. Thus, the 
Hermitean S only contains terms higher than bilinear. 
Now let us assume 



exp(i^)n fc n fc+ i exp(-iS) = ^ 



a, • nini + i 



(Dl) 



Where a\ are constants for fixed k and I. This is im- 
possible, because this expression contains terms at least 
hexa-linear in creation/annihilation operators since S is 
at least quadri-linear. This still holds including £1. 
Next assume 

exp(iS)n fc n fc+ i exp(-iS) oc exp(i/[{n}])c ; t +1 C; (D2) 

and consider (0| c mf exp(iE)n k nk+i exp(— \E)c\ ni |0). 
This had to result in <5 m/ ,i+i<5 mi ,z exp(ig[{n}]). Here, 
/[{n}] is an arbitrary functional of the number oper- 
ators. In general, the functional (/[{n}] differed from 
/[{n}]. However, (0 1 c TO Scj n . |0) = 0, since at least 
two indices of each coefficient £ in S had to coincide 
in order to give a non-zero contribution. But then 
the £ vanish themselves. With the same argument, 
(0|c m S"4. |0) = for all integer n (including zero, 
since mj ^ mf is assumed). Thus, the assumption can- 
not hold. This completes the proof, since including Q in 
the transformation U, the only nonvanishing contribu- 
tions from (0| c TO exp i(Q+E)n k nk+i exp —i(Q+E)c$ ni |0) 
come from terms not including 5 at all. This means, that 
the only contributing terms come from f2 alone. 
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2. Proof that the transform (C6) cannot interchange 
hopping and interaction 

Let us assume, Unjiij+xU^ 1 (in case of XXZ-type 
models) contained a term c\ +1 c k . Applying Eq.(C8), this 
means 



(exp if^^expif^x 

x (exp i£)m f ,i+i(exp if)*^,^!^ 



c ~ c ~ 
mi m f m z 



(D3) 



contained a term ct +1 c fc . H ere, j is fixed and all the m's 
are summed over. Eq.(D3) implies 



m, or 



1. fhi = k, rrif = k + 1 and rrii = rhf 

2. mi = k, mf = k + 1 and fhi — rhf =: m. 
Considering 1.: Eq.(|D3D in this case becomes 



(expif) fc+ ij(expif) 



(D4) 



x (exp i£)m,j+i ( ex P iOfe 4+ 1 ^k ' 
which is proportional to 

(exp -i£)j, m (exp if ) m j+i (exp i^fe+ij (exp '^Otj+i '■ 
x (1 — n m )c k+1 c k . 



Since m is summed over, the 1 in the braces vanishes. 
This is because exp(— if) exp if = 11. Additionally, no 
two-particle hopping must occur. This implies that for 
arbitrary m =/= fh 



(exp if) m j+i (exp if )t j x 

x (expif)fe + ij(expif)t , +1 = 



(D5) 



-Ei 



Vanishing i?i means no occurrence of nearest neighbor 
hopp ing. Hence, E\ is assumed to be nonzero. Since 
(D4)^: =>■ for every fixed j one can find an m so that 
(expif) mj ^ and (expif) mj - + i 7^ 0. Now assume the 
existence of a further m/mso that even (exp if ) iyi j 7^ 0. 
This implied (expif) mj+ i(expif)mj =/= in contradic- 
tion to Eq.( p5| ). Thus, there exists only a single m for 
which (exp if ) m j+i and (exp if ) m ,j don't vanish. But 
this would mean detexpif = 0, which would contradict 
the assumptions. 



Considering 2.: Eq.(D3) now becomes 



(exp if) 



fc+i.j+il 

XC k+l C k C ln C m- 



(exp if )t ,+i (exp -if )j, m (exp if) 



As above, no two-particle hopping must occur. This im- 
plies that for arbitrary m 7^ m 



(exp if ) fe+ i (exp if )£ J+1 



:(expif)m,j(expif)^ , = 0. 



(D6) 



Vanishing E2 means no occurrence of nearest neigh- 
bor hopping. Hence, E2 is assumed to be nonzero. 
Since det exp if 7^ 0, one finds for every fixed j an 
m so that (expif) mj 7^ 0. Now assume a further 
m 7^ m to exist so that (expif)m j 7^ 0. This im plie s 
(expif) mi:) (expif)mj 7^ in contradiction to Eq. flDq ). 
Thus, there exists only a single m for which (exp if ) m j 
is nonzero. This is in contradiction to E2 being nonzero. 

This completes the proof. 
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